arXiv:1501.01795vl [cond-mat.soft] 8 Jan 2015 


Low-density crystals in charged colloids : Comparison with Yukawa theory 


Ioatzin Rios de Anda, 1,2 Antonia Statt, 3,4 Francesco Turci, 1 and C. Patrick Royall 1,5,6 

1 H. H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 1TL, UK 
2 Bristol Centre for Functional Nanomaterials, Centre for Nanoscience 
and Quantum Information, Tyndall Avenue, Bristol BS8 1FD, UK 

3 Institut fur Physik, Johannes Gutenberg-Universitat Mainz, Standinger Weg 9, 55128 Mainz, Germany 

4 Graduate School of Excellence Materials Science in Mainz, Staudinger Weg 9, 55128 Mainz, Germany 

5 Centre for Nanoscience and Quantum Information, Tyndall Avenue, Bristol BS8 1FD, UK 
6 School of Chemistry, University of Bristol, Bristol BS8 ITS, UK 

Charged colloids can behave as Yukawa systems, with similar phase behaviour. Using particle- 
resolved studies, we consider a system with an unusually long Debye screening length which forms 
crystals at low colloid volume fraction <f ~ 0.01. We quantitatively compare this system with the 
Yukawa model and find that its freezing point is compatible with the theoretical prediction but that 
the crystal polymorph is not always that expected. In particular we find body-centred cubic crystals 
where face-centred cubic crystals are expected. 


I. INTRODUCTION 


Among the great successes in describing interactions 
between colloidal particles in suspension is the Derjaguin, 
Landau, Verwey and Overbeek (DLVO) model @6]. 
Apart from short ranged interactions, this model treats 
charged colloids as Yukawa particles. Thus colloids may 
be interpreted within a framework which encompasses a 
great many other systems from the mesons which Yukawa 
was originally interested in m to dusty, or complex plas¬ 
mas mm- Indeed, as shown in Fig. |T] charged col¬ 
loids and complex plasmas both exhibit the same phase 
behaviour : under the right conditions both systems ex¬ 
hibit phases characteristic of the Yukawa system, that is 
to say fluid (F) along with body-centred and face-centred 
cubic crystals (BCC and FCC respectively). 

In the case of colloidal dispersions, there are more com¬ 
ponents than just the colloids : the particles are im¬ 
mersed in a solvent, the electrostatic charge they carry 
is balanced by counter-ions (in colloidal dispersions, one 
typically assumes charge neutrality) and salt ions, not 
to mention the liquid solvent in which the system is 
immersed [48]. One can proceed to an effective one- 
component system where only the colloids are considered 
by integrating out the degrees of freedom of the smaller 
species m- Here the liquid solvent has no impact on 
the equilibrium phase behaviour (it acts to damp the dy¬ 
namics of the colloids, leaving the system as non-inertial 
on most reasonable timescales) |48, . However the effects 
of the ions do need to be integrated out, and this can be 
done using the approach pioneered by Derjagiun, Lan¬ 
dau, Verwey and Overbeek [46]. In addition to the effects 
of the electrostatics, DLVO theory also includes other in¬ 
teractions between the colloids, such as van der Waals 
forces. In our systems, these are not important, because 
the colloids are refractive index matched to the solvent 
which reduces the effects of the van der Waals forces to 
a fraction to the thermal energy VbT. Any residual van 
der Waals effects are suppressed by a polymer layer at 
the surface of the particle. The polymer layer is much 


thinner (< 10/im) than the particle size [51 J and thus 
the short ranged interactions may be treated as a hard 
core [52] . 

In colloidal dispersions, the particle concentration is 
often high enough that steric interactions due to the finite 
particle size can come into play [53], and it is appropriate 
to include the hard core u^r). The hard core Yukawa 
interaction then reads 


u(r) = u hc (r) +u y (r ), 

exp[— n(r — cr)] 


Uy{r) = Cy 


r /c 


(1) 

( 2 ) 


Here, the potential at contact (when the colloids touch) 
is given by 


= 


z 2 Ab 

(1 + kg j 2) 2 a 


( 3 ) 


where Z is the colloid charge, and the inverse Debye 
screening length is given by n=y/4:7r\ Bpiom where pi on 
is the number density of monovalent ions. The Bjerrum 
length 


\ B = /3e 2 /(47re 0 e r ), (4) 

is the distance at which the interaction energy between 
two electronic charges is where e is the electronic 

charge, eo the permittivity of free space, and e r the di¬ 
electric constant. 

In their investigations of hard core Yukawa phase di¬ 
agrams, Hynninen and Dijkstra [53] showed that the 
effect of the hard core was small in the case that freez¬ 
ing occurred at a colloid volume fraction (f> < 0.3, i.e. 
that the colloids are typically far enough apart from one 
another that short-range interactions are irrelevant. Ex¬ 
ample interaction parameters would be (/3e y = 20, ncr = 
1.0). At lower concentrations, therefore, the system can 
be treated as a point Yukawa system, as can also be 
the case with complex plasmas [48] . A point Yukawa 
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FIG. 1. Phase diagram of Yukawa systems. The fluid-solid phase boundary (solid line) is the analytic approximation [Eq. 
[5], the dashed line at large n p denotes that its validity is limited by n p < 10. The approximate position of the BCC-FCC 
crystal boundary is indicated by the dotted line. Narrow coexistence regions are not shown here. Symbols represent various 
crystallization/melting experiments in colloidal dispersions [53] [55] (pink) and complex plasmas [56] (blue). Squares and bullets 
indicate, respectively, BCC and fluid phases (as observed), arrows show the direction in which parameters varied during the 
experiments. Triangles are FCC crystals. Characteristic snapshots of observed fluid and crystalline phases are also shown BE]. 


treatment where the hard core is neglected enables ar¬ 
bitrary Yukawa parameters to be represented in a 2d 
plot [57] EH]. Here we shall use the Yukawa (screened 
Coulomb) coupling parameter T s and scaled inverse De¬ 
bye length k v which are defined below [48] . 

In colloidal dispersions, the dominance of such long- 
ranged interactions means that there is no specific re¬ 
quirement that strongly attractive but short ranged van 
der Waals interactions be suppressed and thus we are not 
limited to sterically stabilised, refractive index matched 
systems. Water and ethanol are popular solvents and sil¬ 
ica and polystyrene are popular materials for the colloids. 
As an aside, this indicates that the melamine particles of¬ 
ten used in complex plasma experiments are in principle 
no different to particles used in colloidal experiments, 
moreover their size ranges overlap [48] . 

While the DLVO theory is only valid in the range that 
linearised Poisson-Boltzmann theory holds (weak electro¬ 
static interactions), higher charging can also be treated 
with a Yukawa interaction by using a renormalised or 
effective charge that is smaller than the physical charge 
on the particles [59] 60 . Thus, providing the effective 
colloid charge can be found, a Yukawa behaviour is recov¬ 
ered. Once the effective charge Z e ^ and Debye length ft -1 
are determined, a number of studies have been made in 
aqueous based systems finding excellent agreement with 
the Yukawa model [6TH64] . 


Behaviour inconsistent with the DLVO model has 
also been seen. The DLVO model recasts interac¬ 
tions between multiple components into a one-component 
Yukawa treatment. Even after accounting for charge 
renormalisation, anomalous behaviour has been ob¬ 
served. In particular the observation of condensation 
like behaviour of the colloids into a “colloidal liquid” 
and “colloidal gas”, with voids appearing in the sys¬ 
tem, was attributed to “like charge attraction” [ 55 ] . 
A variety of explanations have been put forward to ex¬ 
plain this phenomenon, many arriving at the conclusion 
that the effective interaction between the colloids was 
attractive, a surprise for particles with like charge [66L 
[68] . However direct measurement with optical tweezers 
found no evidence of attraction, rather that some ear¬ 
lier measurements (though not the original observation 
by Ito et al. j65]) may have been influenced by arte¬ 
facts maun]. Among the few theoretical explanations 
to have withstood the test of time is that of van Roij 
and coworkers who considered that the entropy of the 
salt ions might drive phase separation to a colloid-rich 
and colloid-poor phase [zn in]. Crucially, in van Roij’s 
treatment, at the two-body colloid-colloid level, a repul¬ 
sion between the particles is maintained : the (repulsive) 
Yukawa form in Eq. [2] is satisfied. The entropic terms 
driving the phase separation do not feature in the (one- 
component) DLVO treatment because degrees of freedom 
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of the small ions are integrated out and captured in a one- 
body term. Quantitative agreement with van Roij’s pre¬ 
dictions was recently found in a phase separating binary 
system whose behaviour would similarly not be expected 
in a pure Yukawa picture [64]. Other deviations from 
DLVO behaviour include ion-colloid decoupling leading 
to a macroscopic electric field which results in extended 
sedimentation profiles mm- 

Of particular interest here is a mechanism originally 
put forward to account for the condensation effects ob¬ 
served by Ito et al. [65b While the two-body term 
between the colloids (Eq. is repulsive, the three body 
interactions induced between three colloids are attractive 
m • Such deviations from two-body behaviour have been 
observed in experiment in the form of a perturbation to 
the fluid structure and indications of non-spherical 
interaction potentials in crystals at higher density m • 
Furthermore simulation work indicates that including the 
three-body terms leads to an increase in the fluid region 
of the phase diagram at the expense of the BCC phase 
[531179]. 

Now the self-dissociation of water means the ionic 
strength is > 10 -7 Mol and additional contributions such 
as counter-ions mean that a Debye screening length of 
ft -1 < 300 nm is typical in experiment. Since parti¬ 
cle resolved studies require colloid sizes of at least a mi¬ 
cron, with aqueous solvents it is difficult to reach condi¬ 
tions where the Debye length is comparable to (or greater 
than) the particle size appropriate for the regime where 
the system behaves purely as a Yukawa system without 
significant contribution from the hard core u^c and other 
e.g. van der Waals short range forces. However in sol¬ 
vents such as cycle hexyl bromide of interest here the 
ionic strength can reach 10 “ 12 Mol so the Debye length 
can be sufficient that micron-sized particles are far apart 
and (point) Yukawa behaviour is found |)4]|55j. In these 
systems, colloidal crystals of exceptionally low volume 
fraction, where the inter particle spacing can run to tens 
of microns have been observed pa uni ihu- These sys¬ 
tems have also been observed to become dynamically ar¬ 
rested, and to fail to crystallise and thus form a glass at 
low colloid density [82j , as indeed have aqueous systems 

m- 


Here we consider the Yukawa parameters associated 
with such “low-density crystals”. Our purpose is to make 
a quantitative comparison between the low-density crys¬ 
tals and the predictions of Yukawa theory in the form 
of the equilibrium phase diagram 0 0 . In partic¬ 
ular we map our experimental data to generic Yukawa 
parameters (T s ,k p ) and compare the state observed in 
experiment to the theoretical prediction. We consider 
the crystal polymorph observed with that predicted. Our 
analysis indicates that while the Yukawa phase diagram 
predicts face-centred cubic crystals for some parameters, 
in our experiments we find body-centred cubic crystals 
only. This contribution is organised as follows : in Sec¬ 
tion [II] we discuss our mapping procedure and assump¬ 
tions, in Section III we outline our experimental tech¬ 


nique. We present our results in Section [IV] and discuss 
these in section El 


II. MAPPING TO YUKAWA THEORY 


The comparison between charged colloids and complex 
plasmas is illustrated in Fig. [l] The main panel in Fig. 
[I] is the Yukawa phase diagram in the T s , k p plane. Here 
we show literature data for complex plasmas [56] and 
colloids [54] [55] . We emphasise that the colloidal parti¬ 
cles illustrated in Fig. [I] are two microns in size. They are 
thus effectively identical to particles used at the smaller 
end of complex plasma experiments [48]. The key dif¬ 
ference is thus the immersing medium, a liquid solvent 
rather than a plasma. 

The freezing line in Yukawa systems is given with rea¬ 
sonable accuracy for k p < 10 [48] by 


T s (/^p) 


106 

1 + fcp + 2 ^ p 


(5) 


where the screened coupling parameter T^ is the 
Yukawa interaction evaluated at the mean inter particle 
separation u y (p~ 3 ) where p is the particle number 
density and k p is a scaled inverse Debye length, given by 
K p =K(j c p [48] [57] [58]. 

Here we assume the colloids take their saturated effec¬ 
tive charge Z^ t . That is to say, the maximum charge 
given under charge renormalisation. An approximate 
value is given by: 


rysat _ ( 2 + K(7 ) cr 

Zeff “ Ab 


(6) 


which represents the effective colloid charge [48]. Thus 
the number density of ions can be estimated as the ef¬ 
fective charge number per colloidal particle due to the 
counter ions and that from salt and background ions p sa i t , 


Pion — ^eff P "b Psalt- (7) 

Now although no salt is added, some background ions 
are present, due for example to solvent self-dissociation. 
Here we treat this contribution p sa i t as a free parame¬ 
ter and determine the scaled screening parameter as out¬ 
lined above. We shall see below that assuming agree¬ 
ment with the Yukawa freezing line enables the value of 
Psait = 10 -10 m ~ 3 which corresponds to 8.3 nMol. Fi¬ 
nally, the evaluation of the Yukawa interaction at the 
mean interparticle separation can be obtained from Eq. 

m 


III. EXPERIMENTAL 

We used poly(methyl methacrylate) (PMMA) colloids 
sterically stabilized with polyhydroxyl steric acid I83ll84l . 






4 


The colloids were labelled with fluorescent rhodamine 
dye to enable fluorescent imaging and had a diameter 
of 2000 nm and polydispersity of 5% as determined with 
static light scattering. Different volume fractions of the 
particles were suspended in cyclohexyl bromide (CHB) 
whose dielectric constant e r = 7.92 and whose refractive 
index is closely matched to that of the colloids enabling 
bulk 3d imaging. A rectangular glass capillary with in¬ 
ner dimensions of 0.10 x 1.00 mm (Vitrocom) was filled 
with the suspensions and sealed on each end with epoxy 
glue. The samples were studied with confocal laser scan¬ 
ning microscopy, CLSM (Leica SP5 fitted with a reso¬ 
nant scanner), with 543 nm excitation using a NA 63x 
oil immersion objective. For qualitative imaging, 2D data 
set was recorded (512 x 512 pixels), whereas for particle 
tracking, a full scan of the capillary in the z direction was 
obtained, providing 3D data sets, where care was taken 
to ensure the pixel size was the same in the three planes. 


IV. RESULTS 


A. Phase behaviour and comparison with Yukawa 
theory 


We show our comparison with Yukawa theory in Fig. 
[2] for the low-density crystals formed in this work. Com¬ 
pared to previous work where the experiments were 
mapped to Yukawa parameters (T s , « p ), we access a new 
region of the phase diagram. We fit the ionic strength 
such that our lowest concentrations 0= 0.0055 and 0.01, 
identified as fluids are consistent with the Yukawa pre¬ 
diction. We note that the path the state points take in 
the (r s , k p ) space is not a straight line or even a smooth 
curve, nor even montonic. This is due to competing ef¬ 
fects. Regadring the ionic strength, the added salt p sa i t 
in Eq. [7] is comparable to or smaller than the counter ion 
contribution Z^p which means that the Debye length 
ft -1 drops as the volume fraction is increased so the 
screening is stronger. However the increase in <p means 
the range at which the screened coupling parameter T s is 
evaluated, p~ 3 5 drops because the particles are (on aver¬ 
age) closer together. The former acts to reduce T s , the 
latter to increase it. We note that the highest three den¬ 
sities (k p ) are predicted to the FCC. We describe how we 
determined the crystal structure in the following section. 

We are now able to quantify the Debye screening 
length ft -1 = 1.9/im and contact potential /3e y = 1110 
around freezing. This corresponds to a effective colloid 
= 850. In particular the Debye length is bigger, 
compared to aqueous systems and indeed to previous 
particle-resolved studies where it was around 0.200-1 p m 

mm- 


B. Identification of crystal structures at low 
packing fractions 

For the analysis of the crystal structures particle track¬ 
ing of the 3d data sets based on [85] was performed. Im¬ 
age manipulation techniques are used in order to enhance 
the contrast and remove the noise from the experimen¬ 
tal datasets. The positions of the particles are identified 
through a maximization of the overlap between seeded lo¬ 
cal gaussian kernels and the intensity profile of the par¬ 
ticles. For this algorithm, an estimation of the size of 
one particle is needed as input parameter, but no initial 
assumptions on the volume fraction are required. The 
resulting coordinates are obtained with a resolution of 1 
pixel (200 nm) in each spatial direction. 

The coordinates extracted via particle tracking are 
used in order to perform an analysis of the crystalline 
phases formed before the sedimentation of the sample. 
For this purpose, we employed the Steinhardt local rota¬ 
tional invariants, also known as bond-orientational order 
parameters. These discriminate between different possi¬ 
ble crystal structures on the basis of spherical harmonics 
and have also been used for the detection of order in com¬ 
plex plasmas [56j[86]. In particular, we consider the lo¬ 
cally averaged order parameters q±,w±, qe for square and 
hexagonal order, where the local average allows to take 
into account the effect of second nearest neighbours and 
to more sharply distinguish between different arrange¬ 
ments (see m for a detailed discussion on the technique). 

For each particle i, we perform a parameter-free detec¬ 
tion of the nearest neighbours via a Voronoi tessellation 
of the sample volume. This provides a list of the Nb(i) 
nearest neighbours over which the local order parameters 
are calculated: 

1 N b (i) 

Qlmii) = Ylm{rij) , (8) 

where Yi m (rij ) are the spherical harmonics. Summing 
over the list of Nb(i) particles identified by the neighbours 
and the particle i itself one obtains the locally averaged 
order parameters 


x N b (i) 

1 (*) = jCr / • \ 

N b( l ) t=0 


qim D J — -jC't / • \ y j (9) 

and summing over all the harmonics we finally get 

Qi(i) = 


47r 


\ 2Z + 1 ^ 

\ rn=—l 


£ \Qlm(i)\ 2 • (10) 


wi(i) = 


E ( mi mo m3 ) 

m 1 +m 2 +m 3 =0 \ " ll " 12 m 3 


i 2 ) 


3/2 


(ii) 
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FIG. 2. Phase diagram of Yukawa system studied in this work in comparison with previous studies [48] . Snapshots in (a) 
(b) (c) correspond to volume fraction 0 of 0.0055, 0.02 and 0.23 respectively and the state points in (r s ,ft p ) representation 
indicated in the main panel. Lines are as in Fig. [I]: thick line is the freezing line [Eq. [5] and the thin dashed line approximately 
describes the BCC-FCC transition. Arrows denote increasing volume fraction. Scale bars = 25 /xm. 


where the term in brackets is the Wigner symbol. 


We focus our analysis on a bulk region discarding top, 
bottom and lateral edges for a thickness of about 2.5 /xm. 
In the lowest density samples, we can discriminate be¬ 
tween a fluid and a solid phase, where the fluid presents 
the characteristics of a layered liquid along the vertical 
z dimension. In Fig. [3] we show the results of the local 
order analysis for a low density system (top row, volume 
fraction </> = 0.015) and a dense sample (bottom row, 
< j) = 0.16). The very limited range of the q 4 order param¬ 
eter allows us to discard the hypothesis of a face centred 
cubic crystal. In order to asses the nature of the solid 
phase, we use an additional order parameter w 4 , particu¬ 
larly suitable for the distinction of hexagonal close packed 
structures (HCP) and FCC from BCC: in Fig. [3jb) we 
show that no peak is detected in the HCP or FCC re¬ 
gion, leading to the identification of the solid phase as a 
BCC phase. We see that the state point at 0 = 0.16 [(c) 
in Fig. [ 2 ] is identified as BCC while the theory predicts 
that it should be FCC. Indeed we analysed all crystalline 
systems and found only BCC crystals. We speculate on 
the causes of this discrepancy below. 


V. DISCUSSION 

Our work extends the range of Yukawa parameters 
(T s , k p ) for particle-resolved studies of colloidal systems. 
We have shown that the formation of low-density crystals 
in these systems [51 ISSUED is compatible with Yukawa 
theory, and that Debye screening lengths can run to many 
microns. The screened coupling parameter is found to 
be T s = 67, the largest value obtained previously was 
T s = 16. We further use our parameterisation to esti¬ 
mate the lowest freezing density attainable in the case 
of no salt (p S ait —> 0). This turns out to be around 
4> = 0.0004 for our system and a Debye length of some 
14 fi m. In principle, therefore it might be possible to 
produce crystals of very much lower density still. 

We now consider possible reasons for the discrepancy 
of the polymorph we identified. One possibility might be 
three-body and higher order effects, which could be sig¬ 
nificant for such long-range interaction potentials. How¬ 
ever as we have noted above simulations predict the op¬ 
posite behaviour, of a suppression of the BCC phase 
relative to the two-body Yukawa system and the for¬ 
mation of FCC instead [53 Eg. Our primary spec¬ 
ulation on the other hand is that the system has not 
yet fully equilibrated. Indeed in some samples, the re¬ 
gion close to the wall appearing more crystalline, indi- 





6 



FIG. 3. (a,d) Three-dimensional rendering of the dense BCC crystal seen from the bottom of the sample for (a) volume fraction 

4> = 0.015 [corresponding to point b in Fig. [ 5 ] and (d) volume fraction 0 = 0.16 [point c in Fig. [ 2 ], with arbitrarily scaled 
particles radii for better visualisation. In panel (d), a close-up of the ordered surface is highlighted, (b-c, e-f) Corresponding 
local bond order parameter diagrams for the bee phase of the two considered volume fractions: in panel (b) the distribution 
P((/ 4 , qo) is centred at moderate values corresponding to a solid region whose nature is determined via the P(q 4 , w±) distribution 
to be of BCC nature. The same procedure is applied in panels (e,f) for the higher volume fraction, showing a higher degree of 
order, particularly in the P{q^^qo) distribution. 


eating a possible crystallisation front beginning at the 
wall of the same cell as has a been observed in the hard 
sphere suspensions [88] . While homogenous crystallisa¬ 
tion for our parameters has yet to be studied in detail, 
we note that in Yukawa systems, although FCC is the 
favoured polymorph, BCC can form first in the Ostwald 
rule of stages and indeed polymorph selection can even 
proceed through the metastable hexagonal close-packed 
polymorph [89]. We believe that the same could oc¬ 
cur here, although we emphasise that the effect of the 
wall will be profound in the crystallisation mechanism 
and thus the homogenous crystallisation studies of Des- 
granges and Delhomme [89] may not hold for our case. 

However we note that our system is confined in a capil¬ 
lary of height 100 /im. Confinement has been considered 
in comparable systems up to four layers [90] . There 
some preference was found for square symmetry which 
might lead to BCC being favoured in our larger system. 
Thus it would be most interesting in the future to in¬ 
vestigate the crystallisation kinetics to see if the system 
indeed showed signs of approaching its FCC bulk equi¬ 
librium state. However we emphasise that these systems 
are stable for around two days m which does place 
some limits on the experimental time window. While no 
change is observed on the timescale of two days (here 


limit our experiments to four hours), after a week the 
system was found to behave as of the contact potential 
/3e y has dropped along with the Debye length. This be¬ 
haviour is inferred from the shift in the freezing bound¬ 
ary to higher colloid volume fraction for older samples 
[Ml [55]. Although desirable, an in-depth study of these 
ageing phenomena has yet to be carried out. Possible 
sources of the change in the system over time include ion 
dissolution from the glass capillary in which the system 
is confined. 
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